Rectification in Y-junctions of Luttinger liquid wires 
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We investigate rectification of a low- frequency ac bias in Y-junctions of one-channel Luttinger 
liquid wires with repulsive electron interaction. Rectification emerges due to three scatterers in 
the wires. We find that it is possible to achieve a higher rectification current in a Y-junction than 
in a single wire with an asymmetric scatterer at the same interaction strength and voltage bias. 
The rectification effect is the strongest in the absence of the time-reversal symmetry. In that case, 
the maximal rectification current can be comparable with the total current ~ e 2 V/h even for low 
voltages, weak scatterers and modest interaction strength. In a certain range of low voltages, the 
rectification current can grow as the voltage decreases. This leads to a bump in the I-V curve. 

PACS numbers: 73.63.Nm,73.40.Ei,71.10.Pm 
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I. INTRODUCTION 

Recently there was much interest in rectification in 
nanoscale systems^ - — Nonlinear mcsoscopic transport 
exhibits much interesting physics such as a peculiar mag- 
netic field dependence of the curren t 16 : 17 and negative 
differential resistance for the rectification current at low 
voltages . 18 ' 19 Another motivation for the investigation of 
mesoscopic diode or ratchet^ 3 - effect comes from possi- 
ble practical applications in nanoelectronics and energy 
conversion. Following the pioneering paper by Christen 
and Biittiker— most attention has focused on a simpler 
case of a Fermi-liquid conductor. At the same time, the 
diode effect requires a combination of spacial asymmetry 
and strong electron interactions in the conductor. Hence, 
one may expect a stronger ratchet current in strongly in- 
teracting Luttinger liquid systems. This expectation has 
been confirmed by a recent study of transport asymme- 
tries in one-channel quantum wires.— ~— Refs . [l8l - l20l have 
focused on a one-channel Luttinger liquid in a linear con- 
ductor in the presence of a single asymmetric scatterer. 
This is the conceptually simplest situation giving rise to 
rectification. At the same time, changing geometry may 
increase asymmetry and hence the rectification current. 
In this paper we consider an asymmetric setup based on a 
Y-junction of three quantum wires with three impurities. 
We show that a stronger diode effect can be achieved in 
such system than in a linear Luttinger liquid and recti- 
fication is possible even in the case of symmetric point 
scatterers. 

We focus on the simplest one-channel Y-junctions. 
More complicated Luttinger liquid junctions, such as Y- 
junctions of single- wall carbon nanotubes, are also of in- 
terest. In particular, it might be easier to make such 
junctions in a reproducible way. 

Y-junctions are among basic elements of electric cir- 
cuits, however, a theoretical investigation of Luttinger 
liquid Y-junctions^— has begun only during the last 
decade. By now, there is a good understanding of linear 
conductance near various fixed points as well as tunnel- 
ing density of states i 33 : 38 In this paper we extend the 
previous research to the problem of transport asymme- 



tries. Specifically, we consider a setup of the type shown 
in Fig. 1. We assume that one of the three terminals is 
kept at zero voltage, ac voltages with amplitudes V and 
jV, 7 ~ 1, are applied to the remaining two terminals. 
In a general dc-current is generated in each of the 

three wires in the junction, Fig. 1. The three currents 
are different but can be computed in a similar way. We 
calculate such rectification dc currents for various types 
of Y-junctions. We focus on the limit of a low-frequency 
ac voltage bias. In order to determine the amplitude of 
the rectification effect in that limit, it is sufficient to find 
the difference of the dc currents at the opposite dc bias 
voltages, i.e., compare the current when the potentials at 
the terminals are time-independent and equal 0, V and 
with the current when the potentials arc 0, — V and 
— jV (cf. Refs. [l8lfl9[ ). This corresponds to a dc cur- 
rent generated by low-frequency square voltage waves. 
We find that in some classes of Y-junctions the rectifica- 
tion current is higher than in a linear wire with the same 
strength of the repulsive electron interaction at the same 
voltage bias. In a certain interval of low voltages the rec- 
tification current exhibits a power-law dependence on the 
bias: / ~ V z . The exponent z can be negative. Clearly, 
such dependence with a negative exponent cannot extend 
all the way to zero voltage as I = at V = 0. Hence, 
the rectification current reaches a maximum at a certain 
voltage. We demonstrate that in junctions without time- 
reversal symmetry, the maximal rectification current can 
be comparable with the total current ~ e 2 V/h even for 
low voltages and modest interaction strength in the wires. 
In particular, such rectification current can be achieved 
in the "island setup", illustrated in Fig. 1, which can be 
experimentally realized in quantum Hall systems^. For 
comparison, in a linear wire, the rectification curren t 18 : 19 
is always much lower than e 2 V/h. The diode effect in Y- 
junctions is not as strong in the presence of time-reversal 
symmetry as in its absence. Still, the maximal rectifica- 
tion current is higher than in a linear conductor. 

The paper is organized as follows: first, after a brief 
qualitative discussion, we describe our setup and for- 
mulate a model. We then review the properties of Y- 
j unctions and derive a general expression for the recti- 
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FIG. 1: Schematic picture of a Y-junction of quantum wires 
with three impurities. Voltage bias ±V is applied to the first 
wire. Wires 2 and 3 are connected to the ground. We calculate 
the dc current in wire 1. One can view AB, CD and EF as 
edges of an integer quantum Hall liquid. This corresponds to 
an "island junction" .— 

fication current. Next, we apply the general formalism 
to the model of three weakly connected wires. Then we 
determine the leading contributions to the rectification 
current at different interaction strengths and matching 
conditions at the junction. We discover several regimes 
with different voltage dependences of the current. Fi- 
nally, we discuss how to build a junction with the maxi- 
mal rectification current. 



II. RECTIFICATION IN MESOSCOPIC 
CONDUCTORS 

We consider a junction of N one-channel wires. In the 
subsequent sections we will specialize to the case of N = 3 
and assume that the system is spin-polarized or, equiv- 
alently, that charge carriers are spinless fermions. In all 
cases we assume that long range Coulomb interaction in 
the conductor is screened by the gates. In such situation 
the applied voltage bias affects not only the current but 
also the charge density in the system. 

Rectification is possible, if a change of the voltage bias 
sign results in a change of not only the sign but also 
absolute value of the current. The latter obviously re- 
quires left-right asymmetry in the system. Asymmetry 
may have two origins: asymmetry caused by the way how 
the bias is applied and geometric asymmetry due to, e.g., 
an asymmetric scatterer . 18 ' 19 In particular, one can con- 
sider a situation in which several terminals are kept at 
zero electric potential and the potential of the other ter- 
minals changes between +V and —V. The charge den- 
sity is different for the opposite voltage signs. This in 
turn affects the current and leads to rectification. Such 
density-driven rectification is possible in a broad range of 
situations including systems without electron interaction 
in the presence of time-reversal symmetry. It has been 
investigated in the context of a carbon nanotube junction 



in the noninteracting-electrons approximation in Ref. UH. 
At a high voltage bias a strong rectification effect was 
found. In this paper, on the other hand, we focus on 
the low- voltage regime and, in particular, the universal 
behavior at low bias. In that limit the density-driven 
rectification effect is of little importance as it results in a 
small rectification current ~ V 2 in a noninteracting sys- 
tem. Besides, the density-driven rectification requires a 
specific way to apply bias. If we are interested, for exam- 
ple, in the transformation of incoming electromagnetic 
radiation into a dc current then clearly the potential os- 
cillates in all terminals. Thus, in this paper we focus on 
the rectification mechanism due to geometric asymmetry. 

Y-junctions have a "built-in" geometric asymmetry. 
Indeed, let us focus on the current I\ in one of the three 
wires connected by the junction. The current I\ = I2 + I3 
equals the sum of the currents in the other two wires. 
Thus, the same current enters the junction through two 
wires and leaves through one only. This means that the 
two sides of the junction are not equivalent since they 
correspond to 1 and 2 wires. Thus, we may expect that 
geometric mechanism of rectification applies to any Y- 
junction. However, as we demonstrate below, the geo- 
metric mechanism only works, if electron interaction is 
present or time-reversal symmetry is broken. 

Indeed, let us consider a time-reversal invariant sys- 
tem of non-interacting spinless fermions in an A-tcrminal 
mesoscopic junction. We calculate the current between 
terminals 1, . . . , K and K + 1, . . . , N. We compare the 
current in the situation when terminals 1, . . . , K are kept 
at the voltage V and terminals K + 1, . . . , N are kept at 
zero bias (case 1) with the current in the situation when 
terminals 1 , . . . , K are kept at zero bias and terminals 
K+ 1, . . . , N are kept at the voltage V (case 2). For non- 
interacting fermions the current reduces to the sum of the 
single particle contributions corresponding to each energy 
in the window < E < V. We thus compare such contri- 
butions for two opposite voltage biases. The wave func- 
tion of an electron incoming from terminal I with the mo- 
mentum k is ipk,i = cxp(-ikxi) + J2m=i Sim cxp(ikx m ), 
where x m > is the coordinate in wire number m 
and Si m is a unitary scattering matrix. The above no- 
tations for the wave function imply that the probabil- 
ity to find an electron with the wave function ipk,i in 
wire number m ^ I is determined by the outgoing wave 
and is proportional to |<Sz m | 2 . For I = m, the proba- 
bility is determined by both the incoming and outgo- 
ing waves. The time-reversal symmetry implies that 

iplj = exp(ikxi) + Yl m =i S*m ex P( — ikx m ) is also a so- 
lution of the Schrodinger equation. The solution is made 
of N incoming waves and one outgoing wave. It can be 
represented as a linear combination of waves ipk,m with 
different m. Hence, the outgoing wave satisfies the equa- 
tion exp(ikxi) = exp(ikxi)^2 m= -LSf m S m i. The absence 
of outgoing waves in the channels with numbers n =/= I 
implies that J2 m S L S mn = 0. Hence, S* = S' 1 = S\ 
i.e., S is a symmetric matrix. Now we can compare the 
contributions to the current from particles with the en- 
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FIG. 2: A schematic picture of a linear wire with two scatter- 
ers of unequal strength. Charge density, averaged over Fiedel 
oscillations, follows a "staircase" profile. The direction of the 
staircase depends on the voltage sign as seen from figures 2a 
and 2b. 



ergy E = h 2 k 2 /2m for two opposite signs of the bias. 
Electrons, incoming from different terminals, are not co- 
herent. Hence, in case 1, h ~ Yli=i 52 m =K+i \Sim\ 2 ■ 

In case 2, I 2 ~ Yd=\ J2 m =K+i \ s ™i\ 2 - The currents arc 
equal from the symmetry of the scattering matrix and 
hence the diode effect is absent. 

In the presence of electron interaction, the rectifica- 
tion effect becomes possible. This can be understood 
already from a simple model with two point scatterers 
of unequal strength in a linear wire (Fig. 2). Similar 
rectification mechanisms operate in more complex junc- 
tions of Luttinger liquids. We will assume that long- 
range Coulomb interactions are screened by the gates and 
thus the charge density depends on the voltage. In Fig. 
2a we consider the situation with the incoming current 
from the left. Backscattering off two impurities results 
in a "staircase" charge density profile (we have averaged 
over Friedel oscillations) . The "staircase" goes up as one 
moves from the right to the left. For the opposite voltage 
sign, the incoming current arrives from the right, Fig. 2b. 
We again have a "staircase" charge density profile but 
now the "staircase" goes down as one moves from the 
right to the left. In the absence of electron interactions, 
the transmission coefficients must be the same in both 
cases and rectification is absent. Let us now consider 
the effect of electron interactions in a simple mean-field 
Hartree picture. Since long-range Coulomb interactions 
are screened, the relation between the charge density p{x) 
and the electric potential W(x) is local. Assuming small 
charge density variations, we thus find W{x) ~ p{x). In- 
coming electrons are scattered by a combined potential of 
the two impurities and the Hartree potential W(x). The 
latter depends on the incoming charge density and hence 
the applied voltage. Hence, for opposite voltage signs, 
electrons feel different backscattering potentials. This 
results in different transmission coefficients and thus rec- 
tification. 

Another mechanism of transport asymmetries comes 
from time-reversal symmetry breaking. The latter is pos- 
sible in the presence of a magnetic field. We illustrate 
asymmetric transport in the absence of the time-reversal 
symmetry with a model system depicted in Fig. 3. While 
it is not a realistic model of any mesoscopic conductor, it 
is very simple and exhibits strong transport asymmetries. 
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FIG. 3: Charge transport through the box exhibits no asym- 
metries at zero magnetic field, Figs. 3a,b, and is asymmetric 
at a finite field in shaded areas, Figs. 3c, d. 



We assume no electron interactions in the model. 

We consider a box with four holes in its sides. Each 
hole corresponds to a terminal in a more realistic de- 
scription of a junction. We assume that all charge carri- 
ers, entering the box through the holes, have exactly the 
same speed, perpendicular to the wall in which the hole 
is made. The speed v = \/—2meV, where V plays a role, 
similar to the electrostatic potential of a terminal in a 
realistic junction. At V = the speed v — and hence 
no current is injected. In a Y-junction, charge density, 
injected from each terminal, depends on its electrostatic 
potential. A comparison between opposite voltage signs 
in a realistic junction corresponds in our model to the 
comparison between the situations with the charge in- 
jected through the right and left holes (Figs. 3a, c) and 
through the upper and lower holes (Figs. 3b, d). In the 
absence of a magnetic field, the current through each 
hole is zero in both situations: currents of the particles, 
injected through the opposite holes, cancel as shown in 
Figs. 3a, b. Let us now turn on a magnetic field in shaded 
areas (Figs. 3c, d). If the charge carriers are injected from 
the upper and lower holes then the current remains zero 
even in the presence of a magnetic field (Fig. 3d). Indeed, 
the carriers never enter the region with the field. How- 
ever, if the carriers arrive from the left and right holes 
then they are deflected by the field (Fig. 3c) and the elec- 
tric currents are nonzero in all four holes. This illustrates 
how magnetic field can result in transport asymmetry. 

The above arguments make it plausible that the rec- 
tification effect is the strongest if both strong electron 
interaction is present and time-reversal symmetry is bro- 
ken. This is confirmed by our calculations below for 
Y-junctions of quantum wires with and without time- 
reversal symmetry. 
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III. Y-JUNCTIONS 

In this section we formulate our model and review basic 
properties of Y-junctions. 

We consider a Y-junction with the action 

fc=l,2,3 fe=l,2,3 

where are the actions of three uniform wires and Tk 
describe three impurities in the wires close to the junc- 
tion point x = 0. The above action must be supplied by 
matching conditions for three wires. They will be dis- 
cussed below. The action of a uniform wire is given by 
the equation 

poo 

L k = / dx\ti$ (d t - v F d x )i> k j + (d t + v F d x )^ k 
Jo 

-v^f^+^^of], (2) 

where ipi o are the operators of incoming and outgoing 
chiral electron fields, Vp is the effective Fermi velocity, 
vp\ defines the interaction strength. We set h = 1 and 
the electron charge e = 1 in most of the following text. 
We assume that the long range part of the Coulomb force 
is screened by the gates. Hence, the electron density 
depends on the voltage bias. The impurity Hamiltonians 

/>oo 

T k = dxU k {x)[i>fijl + h.c], (3) 
Jo 

where Uk{x) is the potential of the impurity in wire num- 
ber k, the impurity being located close to the junction at 
x = 0. We only take the backscattering part of the impu- 
rity Hamiltonian into account as the forward scattering 
terms do not affect our results. We assume the same 
interaction strength and Fermi velocity in each wire. 

In order to treat the case of strong interaction it is 
convenient to bosonize^ the action in terms of the chiral 
fields (j)Q j such that V'o/i = ^o/i e x p(±*&F£ + i(/>Q/j), 
where the opposite signs should be chosen for the in- and 
out-fields, the commutator 

[<l>o/i(y), d x<f>o/i( x )] = T2TTi5(y - x)6 u (4) 

and k F plays the role of the effective Fermi momentum 
and determines the average charge density in the wires. 
Fj^o are Klein factors, necessary to ensure the proper 
Fermi commutation relations. The local densities of in- 
and out-moving particles in point x are Po/i = (k F ± 
d x cj)o/i)/2iT. The action now assumes the form 

L k = J dx±-[d x <t> k (d t - v F d x )$ 

+ d x ^ k (-d t - v F d x )$> - —(Mi " dx4>of\\ (5) 



T k = £ U n (F^Fj) n exp(in[<ft(x - $)-<p k {x = 0)})+h.c, 

n 

where U n = U n cxp(ia n ), with real U n and a n , are of the 
order of the Fourier components of the asymmetric po- 
tential, kp J ex.\>{i2nkpx)U {x)dx. Note that a n can be 
nonzero even for a symmetric potential U(x) in contrast 
to the situation considered in Refs. [l8lfl9l For example, 
for U ~ 5{x— xq), Cii = 2kpXQ. The above expression for 
the backscattering operators Tk includes multi-particle 
backscattering processes. 18 ' 19 Such multi-particle contri- 
butions to the action are inevitably generated under the 
action of the renormalization group by the interplay of a 
short-range Coulomb interaction and impurity potential 
(see Ref. [l9| for a discussion) . All prefactors U n have the 
same order of magnitude and are proportional to U(x). 
We assume that backscattering amplitudes U n have the 
same order of magnitude in all three wires. The use of 
the fields <fi at x — in the backscattering operators is 
justified, if the distance from the impurities to the junc- 
tion is lower than the scale Hvf / {eV) set by the voltage 
bias. We assume that U n are sufficiently small so that a 
perturbative expansion in powers of U n can be developed 
for the calculation of the current. The conditions on U n 
will be formulated below. 

Following the standard notation conventions, the 
Hamiltonian, corresponding to Eq. ([5]), can be written 
as 

Hk = l^j ^[5(^^i+^^) 2 +i(^0f-5 x ^) 2 ], (7) 

whe re the d imcnsionlcss interaction strength g = 
+ 2X/tt < 1 and v = v F /g. 
The above action alone is not enough to describe the 
system and matching conditions at the junction are nec- 
essary. As we will see, matching conditions are sub- 
ject to several restrictions. The most general match- 
ing condition has the form F(<p k (x = 0), d x (f) k Q (x = 
0), d x (f>j , , . . ) = 0, where F is an arbitrary function. 
Below we will focus on the situation in which the system 
is close to a fixed point. Then only most relevant opera- 
tors should be kept in the matching conditions and hence 
all derivatives of the Bose-fields <fi k can be neglected. 
Next, we note that the shift of any field <j>io by a con- 
stant <fi(x) — > 4>(x) + C is a gauge transformation that 
does not change the physics of the system. Hence, the 
matching conditions must be invariant with respect to 
such gauge transformations. This singles out the bound- 
ary conditions of the forrn^ 

<p k o=Y,M kj 4, (8) 

3 

where M is a matrix with real matrix elements. With 
this matching condition, it is easy to solve the equations 
of motion for the chiral fields tpQ in the absence of the in- 
teraction (A = 0) and impurities, for arbitrary initial con- 
ditions for fields <fq Q . Substituting the solution into the 
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commutation relations (J4j) for in- and outgoing fields, one 
finds that they are compatible, provided that the matrix 
M is orthogonal. The same condition for the matching 
matrix must hold for a general case with an arbitrary in- 
teraction strength— since the problem can be reduced to 
a model of non-interacting bosons by diagonalizing Eq. 
(|7|) in each wire (see the next section). 

Thus, we established that M is a real orthogonal ma- 
trix. We next demonstrate that the sum of the ele- 
ments of each of its rows and columns equals one— 
It is again convenient to consider the situation with- 
out electron interaction and tunneling. In that case the 
currents of in- and outgoing electrons in each wire are 
Ij, = ^fVFpj, = ±dt<j>^, /2n. Charge conservation 
implies that the sum of the currents in all wires is zero 
at x = 0. Taking into account Eq. <JSj> one finds that 
J2k M kj = !• Multiplying Eq. © by M -1 = M T 
and repeating the same argument one also finds that 
Mkj = 1. Again, electron interactions do not affect 
the above result—. 

The origin of the matching conditions can be under- 
stood if one notes that they can be imposed by adding 
to the Hamiltonian a term of the form 



A]Tcos(r#* (x = 0) Y^M, y y r .r = 0)]), (9) 



where A is a large constant. Each cosine tends to keep 
its argument at zero. This can be achieved simultane- 
ously for each cosine only if their arguments commute. 
This happens for a unitary M. The cosines describe tun- 
neling between different wires. Thus, charge conserva- 
tion implies Y]^ Mkj = 1 • This provides an alternative 
derivation of the matching conditions. 

Since we plan to separately investigate time-reversal 
invariant and non-invariant systems, we next need to 
determine what matching conditions satisfy the time- 
reversal symmetry. This is easy as the time-reversal 
transformation corresponds to the change of variables 
4>i/o Wo 1 1- From Eq- ® we then see that in time- 
reversal- invariant systems M = M . Taking into ac- 
count that M is orthogonal we conclude that in time- 
reversal invariant systems it is also symmetric— The ac- 
tion (fTJ) is always time-reversal invariant. Thus, the be- 
havior of the whole junction with respect to time reversal 
is determined solely by the symmetry of the matching 
matrix M. 

At this point we are in the position to give a full clas- 
sification of fixed-point matching matrices.— For time- 
reversal invariant systems, it is convenient to use a 
parametrization from Ref. lU We discover three pos- 
sibilities with time reversal symmetry: 



M = 




M 




M 




(11) 



(12) 



where b = a(a + 1)/(1 + a + a 2 ), a = (a + 1)/(1 + a + 
a 2 ) and c = — a/(l + a + a 2 ). Eq. (fTO]) corresponds 
to three disconnected wires. The cases of a = 0, — 1, oo 
correspond to a junction of two wires and a detached 
third wire. The physics is the same as for one linear wire 
and will not be discussed below. 

In the absence of time-reversal symmetry we use the 
parametrization from Ref. \3§t 



M 




(13) 



(10) 



where a = (1 + 2cos0)/3, 6(c) = (1 - cos 9 ± \/3 sin 6) /3 
and 7^ 0,7r. The case of 6 = 1, a = c = corresponds 
to the island setup, depicted in Fig. 1. Note that the 
matching matrix (|11| has the form (fT3|) with 9 = it so we 
will consider it with the case of no time-reversal symme- 
try. 

In general, matching matrices contain negative ma- 
trix elements. This means that the incoming current in 
one wire may suppress outgoing current in that or other 
wires. Such situation is possible in the presence of An- 
dreev scattering, if, for example, a part of the junction 
is superconducting. Negative matrix elements were also 
predicted in Y-shaped beam splitters for cold bosonic 
atoms— Note that for noninteracting electrons no neg- 
ative matrix elements are allowec— since in the absence 
of interaction the elements of the matrix M reduce to 
the squares of the absolute values of the elements of the 
scattering matrix (cf. Sec. II). 



IV. ELECTRIC CURRENT 

In order to calculate the current we need to include 
a voltage bias. We model Fermi-liquid leads by assum- 
ing that electron interaction in the wires is zero at large 
distances x from the junction— ~— We will assume that 
leads 2 and 3 are kept at a zero voltage. The bias volt- 
age ±V is applied to lead 1. The results do not change 
in a more general model with the potential ±V applied 
to lead 1, the potential ±jV, 7 ~ 1, applied to lead 2, 
and a zero potential applied to lead 3. We will use the 
language of zero 7 since it is simpler. Our approach can 
be easily extended to a general 7. 

The calculations will be based on a renormalization 
group approach with the voltage V playing the role of 
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the infrared cutoff. Thus, we will assume that the tem- 
perature T < V . We would like to emphasize that the 
role of the voltage docs not reduce solely to that of a 
cut-off; otherwise the density-driven rectification would 
be lost (see a discussion after Eq. (fT8|) ). It is well known 
that in Luttinger liquids a voltage bias can play a more 
prominent role than just a cut-off (see, e.g., Refs. [50l — 
1531 ). At the same time, the leading contribution to the 
total current at zero temperature can be estimated by 
simply setting the renormalization group cutoff to the 
value of the voltage^. Below we find that the rectifi- 
cation current is comparable with the total current in 
certain regimes in the absence of the time-reversal sym- 
metry. This certainly means that in those regimes the 
rectification current can be computed by assuming that 
the voltage plays the role of a cutoff only (we however do 
not make such an assumption). 

For a general 7, we do not expect the current to de- 
pend significantly on the temperature at T < V. Such 
dependence can emerge at particular values of 7. In par- 
ticular, in the island setup (Fig. I), this happens at 
7 = 0,1. This can be seen from the Keldysh perturba- 
tion theory^*. The infrared cutoff in some integrals in 
the perturbation expansion is set by max(T, 7V) and in 
some others by max(T, [1 — 7] V). Hence, at 7 = 0, 1 our 
approach applies only for V ~ T in the island setup. 

We want to find the current in the first wire. At low 
frequencies the current conserves and hence it is sufficient 
to find the current in lead 1. It is given by the sum 
of the chiral incoming and outgoing currents ±.vfPi/o- 
The incoming current can be found from the Landaucr 
formula and is linear in voltage. It does not contribute to 
the rectification current. We thus focus on the outgoing 
current in lead 1. It can be found with a generalization 
of the approach of Ref. Our approach is also related 
to that of Ref. US 

In what follows we set the temperature T = to sim- 
plify notations. Our method can be easily generalized to 
finite temperatures. 

Let us introduce an auxiliary field <f>\ (x) : 

4>x{x) =<j) x {x),x> 0;</>i(.t) = ^M lk 4> k I {-x),x < 0. 

fc 

(14) 

The field 4>\ satisfies a simple matching condition 
</>i(+0) = (pi (— 0). Thus, the auxiliary field is a chi- 
ral field propagating through the junction. In the sta- 
tionary regime the average time-derivative of any op- 
erator is zero. Let us now consider the operator O = 
f_ a dxd x <fii/2ir J where the integration extends between 
points taken in the nonintcracting leads. Its mean- 
ing is the charge carried by the mode 4>\. The equa- 
tion = (dO/dt) = i([H,0]), where the Hamiltonian 
H = J2k(Hk + ?fc) describes the whole system including 
the leads, wires and scatterers, reduces to the following 
relation: 

v F {p l {a)) = (vp^MikP^a)) +i([J2Tk,0]). (15) 

fc fc 



The left hand side is the outgoing current we want to 
find. The first term in the right hand side is linear in 
the voltage bias and cannot contribute to the rectifica- 
tion current. Thus, we have to calculate only the sec- 
ond contribution to the right hand side. In other words, 
our problem reduces to the calculation of the average 
backscattering current whose operator equals / = h > 
where 

h =i[Ui(F$F?) n exp(ia; + in{4> k M - 4(0)}) + h.c, 6} 

=i\Ui(F^Ff) n exp(iai + in^A/" 1 )^ (0) - </£(0)}) 

p 

+ h.c.,6] 

^nUtiFgFfyexpiiat+in^M-^kpifeW - ^(0)]) 

p 

x [(M- 1 )*! - S kl ] + h.c. . (16) 

Note that for any I the current operator expresses via all 
backscattering operators in all three wires. 

In order to find the average current (/;), we will as- 
sume that the backscattering operators T k are absent at 
t = — 00 and are then gradually turned on. Without the 
operators T k the system can be viewed as an equilibrium 
one in the ground state of an appropriate effective Hamil- 
tonian. In order to find it we introduce another auxiliary 
field with the structure, similar to <j)\: 

4>i{x) =ct> 1 I (-x),x<O;0 1 (x) = Y^{M- 1 )ik4>o{x),x>Q. 

fc 

(17) 

In the absence of the backscattering operators T k , the op- 
erator A = dxd x (f>i 1 2ir commutes with the Hamilto- 
nian. It can be understood as the charge of a chiral mode 
propagating through the junction. In other words, it is 
an additive integral of motion. Hence, the system can be 
described by a Gibbs distribution with an appropriate 
thermodynamic potential conjugated to A. The physical 
meaning of that thermodynamic potential is the applied 
voltage bias V. At zero temperature, one finds that the 
system is in the ground state of the effective Hamiltonian 
W = H — VA, where H is the actual Hamiltonian of the 
Y-junction. 

The current, i.e., the average of the sum of the oper- 
ators i[Tk,(J], can be now calculated with the Keldysh 
technique^. It is convenient to apply the interaction 
representation H — > H — VA. The interaction represen- 
tation introduces time-dependence into the operators Tk 
and according to the rule: 

exp(*n[^(A/- 1 ) fcp ^(0) - ^(0)]) -> 
p 

exp(iVAt) exp(m[c^(0) - ^ M fcp $(0)]) exp(-iVii) 

p 

= exp( l n[J](A/- 1 ) fcp ^(0) - $,(0) + Vt(M kl - 5 kl )}). 
p 

(18) 
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After such time-dependence is added into all backscat- 
tcring operators, the contribution of the form —V A can 
be removed from the action by a linear shift of all fields 
4>qij- This does not mean, however, that the voltage 
would only enter the action through the time-dependence 
of Tfe. Indeed, since the charge density depends on the 
voltage, the amplitudes Uk of the operators may get 
corrections proportional to the small voltage bias. This 
effect is discussed in Ref. Ojl 

The current can now be found with a perturbative ex- 
pansion in powers of Uk from the standard expression 



h = (0|S(-oc, 0)7,5(0, -oo)|0), 



(19) 



where S(t,t') is the evolution operator from time t' to 
time t and |0) is the ground state of the effective Hamil- 
tonian H — V A. To complete the calculation, we need 
the Green functions determined by the quadratic part 
of the action <j2j> and the matching conditions <JSj> . They 
have been found in Ref. |38|. A Bogoliubov transformation 
of the form cf> /i = (1/2^/77) [(1 + g)tj> /i + (1 - 9)h/o] 
in each wire allows one to obtain free chiral fields 4>i/o 
with correct Luttinger liquid commutation relations^ 
The correlation function of the incoming fields at x = 
is given by the free particle relation 

(#(f)#(tO) - -hv In(i[* - t']/rc + S), (20) 

where t c is the ultraviolet cutoff time of the order of 
the inverse bandwidth, 5 is infinitesimal. The matching 
conditions for the new fields have the forrn^ 



l k {x 



0) 



M kp $(x = 0): 



(21) 



M = [(1 + g) - (1 - g)M]- 1 [(l + g)M - (1 - g)]. (22) 

In the time-reversal invariant case M = M. In the ab- 
sence of the time-reversal symmetry, M has the same 
general structure (fT5|) as M and satisfies the same set of 
constraints but the matrix elements arc different: 



3g 2 - 1 + (3 5 2 + l)cosfl _ 
3(1 +g 2 + (g 2 - 1) cos 6>)' 



6(c) = 



2(1 - cos 6> ± V3gsin6 
3(l+3 2 + ( 5 2 -l)cos( 



(23) 



(24) 



At this point one can write an expression for the cur- 
rent as an expansion in powers of backscattering ampli- 
tudes {/;. Evaluation of the terms of that expansion is 
technically difficult and not very informative. Indeed, 
neither the amplitudes nor the ultra-violet cutoff are 
known exactly. As a result, it is only possible to estimate 
the order of magnitude of each contribution. At the same 
time, such estimation can be performed without explic- 
itly calculating integrals over the Keldysh contour and 
will be sufficient to find the leading power dependence 
of the rectification current on the bias voltage at small 
Ui and low voltages near each of the fixed points. Such 
estimation will be the focus of the remaining sections. 



V. SCALING DIMENSIONS 

It is convenient to use a renormalization group point of 
view for the calculation of the current.— In the renormal- 
ization group procedure, the coefficient Ui in the Hamil- 
tonian and the operators Ii scale as E* 1 ' 1 , where E is 
the energy scale and zi the scaling dimension. At the 
energy scale E ~ V the renormalization group proce- 
dure stops. Different perturbative contributions to the 
current can be estimated from the scattering theory 18 i 19 
The current can be expressed as an infinite sum of the 
contributions, proportional to the products of different 
combinations of Uf. 



I 



const x 



(25) 



Strictly speaking, it is not enough to include only 
backscattering operators from the bare tunneling Hamil- 
tonian ([6]) in the expansion (|25|). All operators gener- 
ated by the renormalization group must also be included. 
Their general structure is 

Wi =Wi exp{i ai + i n k {(b k o(x = 0) - <fi(x = 0)]) 

k 

=Wi exp{i ai + iy/gJ2 J2 n kWk P ~ &kp]<fc)> (26) 

k p 

For simplicity we use the same notation ai for the phases 
of Wi and Uk- Note that it is always sufficient to keep 
only two different values of k in the sum in Eq. (|26|) . since 
~5Zk=i\Mkp — Skp] = 0. This means that each possible 
operator W p can be generated from a product of two 
operators T). ([6]) and hence W p ~ U 2 . Since Klein factors 
do not change the scaling dimensions of Wi , wc omit them 
in Eq. ([26]) . 

To complete the calculation of the current we need to 
find zi. A straightforward calculation based on Eq. (f2"0)) 
yields for the operator (|26|) with an arbitrary choice of 

Uk- 



2/ 



(27) 



In the absence of the time-reversal symmetry the above 
expression greatly simplifies: 

zi = g(l — a)[n\ + n 2 + n 2 — nxn% — nin% — n^n\\. (28) 

In the low- voltage regime, the main contribution comes 
from the operators with the smallest zi. One easily sees 
that it equals 



Zmin = g(l-a) 



(29) 



and is achieved, if two of the coefficients rij = with the 
remaining one being ±1. 

In the time-reversal invariant system, the expression 
for the scaling dimensions is more complicated. It can 
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be simplified by setting m = n\ — ri2 and k = ni 
Using the matching matrix (|12p . one finds: 



n 3 . 



z = g- 



[k(b + c) +mc] 
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1 + a + a 2 



(ka-m) 2 . (30) 



We will assume for simplicity that a is irrational. We 
will thus avoid situations with z — at some k and m. 
Certainly, such situations correspond to the operators 
Wi with no ^-dependence. Such operators cannot affect 
transport. 

The renormalization group procedure only applies if all 
Ui are small at the initial energy scale ~ 1 /r c and remain 
small up to the scale V . In the case without time-reversal 
symmetry, this condition can be easily expressed in terms 

of — 1111 n • 



UiT c < const (Vr c y 



constV 1_ff ( 1_8) . 



(31) 



The above equation assumes that 1 — z min > 0, i.e., Ui 
is relevant. If 1 — z min < then the only restriction on 
Ui is Uit c < 1. What happens in the presence of the 
time-reversal symmetry will be addressed in subsequent 
sections. 



VI. RECTIFICATION CURRENT 

The above discussion applies to all contributions to 
the current, both even and odd in the voltage bias. We 
are interested in the even contribution, i.e., the rectifica- 
tion current I r (V) = [I(V) + 7(-V)]/2. What terms of 
the perturbation expansion contribute to the rectification 
current? The answer to this question can be obtained 
from symmetry considerations. 

In what follows we will need the commutation relations 
for the Klein factors in . Using the commutation rela- 

tions V#(0)$(e) = -V?(e)V#(0) and [<fi(x),<ft(y)\ = 
— wrsign(;r — y), an expression for ipQ in terms of 0j and 
the Baker-Hausdorff formula, one finds 



F$F* = exp(*7r[l - M kk ])Fj F ( 



k pfet 

o ■ 



(32) 



Next, some terms of the perturbative expansion (|25p 
arc zero identically as only certain combinations of the 
vertex operators in T k and // produce non-zero results 
after averaging with respect to the quadratic part of the 
action ([5]). The condition is well known: 



(JJ cxp(i ^ c ifc 0j(O))) ^ only if for each k ^ c ik 



k 



I 



0. 



(33) 

Let us now apply the above results to possible con- 
tributions of different orders of the perturbation theory 
to the rectification current. Strong limitations emerge 
for second order contributions. Indeed, Eq. (|33|) implies 
that any second order contribution, proportional to an 
operator Ui, must also contain its Hermitian conjugated 
operator W . Hence, the phase factors exp(±iaj) can- 
cel each other and drop out from the expression for the 



current. Thus, we can just set ai = 0. Let us now 
compare the currents at the bias voltages V and —V. 
In the case with the bias —V, we make the change of 
variables 4>i/o ~^ —( ri/o m This transformation does not 
affect the form of the quadratic part of the action ([5]) 
and the linear matching conditions. Let us also change 
the order of the Klein factors in the tunneling and cur- 
rent operators using the commutation relation (|32p and 
redefine F — » F+, — » F. Since Ui pairs up with l/j in 
the nonzero second-order contribution to the current, one 
can easily see that the phase factors exp(±z7r[l — M kk \) 
in the commutation relation (|32[) drop out. At the 
same time, if we omit the phase factors exp(±«a;) and 
cxp(±i7r[l — Mkk]) from the action and current opera- 
tors, we discover that the action assumes precisely the 
same form as for the voltage bias V. On the other hand, 
the current operator // (|16p changes its sign. This means 
that the average of Ii also changes its sign and hence 
the second order term of the order UfV 22,1-1 does not 
contribute to the rectification current. 

The above argument assumes that the backscattering 
amplitudes Ui are independent of the voltage bias. Since 
the injected charge density depends on the voltage in 
our setup, Ui can exhibit a weak linear dependence on 
the bias voltage. This means that second order terms 
can in fact lead to a "density-driven" rectification effect 
but its amplitude is suppressed by an additional factor 
Vt c , where the ultra-violet cutoff r c is of the order of the 
inverse band width. Thus, the second order contribution 
to the ratchet current scales as 



I' 2 



UfV 2zi . 



(34) 



Beyond the second order, no additional general restric- 
tions on possible contributions to the rectification current 
can be derived. Particular contributions to the ratchet 
current can disappear at particular values of the phases 
ai and interaction strength g. 



VII. THREE WEAKLY CONNECTED WIRES 

So far we discussed general features of Y-junctions. We 
now determine the rectification current in different se- 
tups. As a warm-up exercise we investigate the simplest 
situation of three almost disconnected wires, Fig. 4. This 
situation corresponds to the matching matrix (jlOp . The 
rectification current is much weaker in this limit than for 
other matching conditions. At the same time, the quali- 
tative picture is quite similar. 

There is no current at all in our original model with 
the action (JT|) and backscattering operators (J3J) as it de- 
scribes three disconnected conductors. We thus modify 
the model in this section. Instead of the backscattering 
operators ^ we consider weak tunneling between the 
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FIG. 4: Schematic picture of a Y-junction of three weakly 
connected wires. Voltage bias ±V is applied to wire 1, and 
bias ±yV to wire 2. Wire 3 is connected to the ground. We 
calculate the dc current in wire 1. 

wires. The tunneling operators have the form 

3 

fe=i 

T k =Fl +1 F k U k e* a * exp(z0 fc (O) - i&+i(0)) + h.c, (35) 

where </>fc(0) = 0/(0) = </>q(0), ^fc are Klein factors and 
we use the convention 3 + 1 = 1. The commutation 
relations for the Klein factors are F k F q = —F g Fk, k ^ q. 
We assume that wires 3 is kept at zero voltage. The 
chemical potential of the first wire is ±V . The second 
wire experiences the bias ±"/V. Note a different meaning 
of the amplitudes U k and phases a k from the model with 
three weak scatterers discussed above. In Eq. (|35|) . the 
phases a k describe the Aharonov-Bohm effect due to the 
magnetic field through the junction. On the other hand, 
in the model with weak scatterers, all information about 
the magnetic flux is contained in the matching conditions. 
The voltage bias can be included in the action in the 
form of the time-dependence of the tunneling operators, 
similar to the discussion above. We will not need explicit 
expressions below. The operator of the current, tunneling 
between wire 1 and wires 2 and 3, is 

I = iT x - %T\ - iT 3 + iTl . (36) 

The current can be estimated from a renormalization 
group procedure that stops at the scale E ~ V. At that 
scale Uk — > UkiVTc) 1 ' 9 ^ 1 , where g < 1 characterizes 
the interaction strength. The second order contribution 
to the rectification current can now be easily found and 
scales as 

jl2) „ u 2 v */g 7 2/ g +l i (37) 

where U ~ Ui < 1/t c . 

The n-th order contribution to the current cannot ex- 
ceed - (UT c ) n (VT c ) n / 9 - n+1 /T c , n > 3. Thus, the leading 
higher-order contribution corresponds to n = 3, 

/(3) „ ^3^3/9-2^3/^ (3g) 



Interestingly, at low voltages and g > 1/2, the third order 
contribution exceeds the second order contribution. 

The above discussion ignored the issue of the time- 
reversal symmetry. Depending on the presence or ab- 
sence of the magnetic field through the junction, the sys- 
tem can have or have no time-reversal symmetry 28 ! 29 It 
is instructive to investigate the effect of the symmetry 
breaking. In the presence of the time-reversal symme- 
try the magnetic flux is zero and all a k = 0. Indeed, 
the time-reversal transformation can be represented as 
4> k j -> 0o,0o ~> <t> k n a k -> -otk- Since <$(0) = 0&(O), 
we conclude that a k = — ot k = at zero magnetic field. 
On the other hand, in the presence of the magnetic field 
a k £ 0. 

The second order contribution to the rectification cur- 
rent does not depend on the phases otk- Let us com- 
pare the third order contributions for different values 
of a k . Nonzero third order contributions to the cur- 
rent originate only from the product of all three op- 
erators T k or all three operators in the perturba- 
tive expansion of the Keldysh expression for the cur- 
rent (fl~9|) . The first contribution, l[ 3 \ is proportional 

to exp(iJ2 k ct k ), i] (V) = exp(iJ2k a k)Ji, and the 
second contribution is proportional to exp(— i^2 k ct k ), 

I^(y) = exp(—iJ2 k a k)J2- Let us change the voltage 
sign. Simultaneously, we change 4>o/i ~^ — ^o/i' cnan g e 
the order of Klein factors in all operators and redefine 
F k — > Fu,f1 — > F k . This transformation is equivalent to 
changing the sign of a k in all operators and simultane- 
ously changing the overall sign of the tunneling Hamil- 
tonian (|35l) but not the sign of the current operator. 

Hence, the third order current l[ S) (-V) + 4 3) (-V) = 
cxp(—iY^ k a k)Ji + cxp(i J2k a k) J%- We see that the 
current is an even function of the voltage in the pres- 
ence of the time-reversal symmetry, i.e., at ^2 k a k = 0. 
The third order current is an odd function of the volt- 
age and does not contribute to the rectification effect, if 
Sfc a fc = i 7r /2. Thus, for three weakly connected wires, 
magnetic field suppresses rectification. We will see below 
that typically it has an opposite effect and rectification 
is stronger in the absence of the time-reversal symme- 
try. Note that in the above example, the ratchet current 
returns to its maximal value at y , a k = tt. 

Let us finally compare our results with the case of a 
linear wire with a high asymmetric potential barrier or, 
equivalently, two weakly connected wires. In that case 
the rectification current is given^ by Eq. (|37|) . Thus, at 
g > 1/2 the rectification effect is stronger in a Y-junction 
than in a linear wire with the same interaction strength. 



VIII. RECTIFICATION EFFECT IN THE 
ABSENCE OF THE TIME-REVERSAL 
SYMMETRY 

Here we address matching matrices of the form f) 13|) . 
In particular, the results of this section are relevant for 
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the island setup, Fig. 1. Our analysis also applies to the 
time-reversal-invariant situation with the matching ma- 
trix (fTTj) . Below we find that at sufficiently low voltages 
the rectification current scales as Eq. ([55]) with U being 
the impurity potential strength, V the voltage, g < 1 the 
Luttingcr liquid parameter, characterizing electrostatic 
interactions, r c the ultraviolet cutoff of the order of the 
inverse band width, and a is determined by the matching 
conditions according to Eq. (|23|) . 

Similar to the previous section we need to compare 
the leading second and third order contributions to the 
rectification current. They are determined by the opera- 
tors of the form Tf, ~ exp(i[<f)j — 4>o\) with the minimal 
scaling dimension (|29l) . The leading second order con- 
tribution scales as ~ (t/r c ) 2 (VT c ) 29 ( 1_a )/r c , where 
5, is given by Eq. (f2"3"| and satisfies —1/3 < a < 1. The 
leading third order contribution comes from the terms in 
the perturbation expansion proportional to the product 
of all three Tfc (or all three T^). It scales as 



low voltages and achieves its maximal value of the order 
of e 2 V/h at V = V* ~ tfV[i-flU-5)]. i n the last equa- 
tion we omitted a power of r c as this does not lead to a 
confusion. 

Finally, let us compare the result with the case of 
a linear wire with the same interaction strength g in 
the presence of an asymmetric impurity. In that case 
the maximal rectification current, calculated in the per- 
turbative regime in Ref. [3 scales as V 1+3g < V at 
g < 1/3 (we omit powers of r c ). The rectification current 
does not exceed V 2 /t c at g > 1/3. In a Y-junction, at 
1— 0(1 — a) > the rectification current ~ V at V ~ V* is 
always greater than in a linear wire. On the other hand, 
a negative 1 — g(l — a) implies g > 3/4. The maximal 
rectification current corresponds then to U ~ l/r c and 
scales as ~ V 39 ^~ a ' 2 > y 2 . Thus, the maximal ratio 
of the rectification and total currents is always higher in 
a Y-junction than in a linear wire in the Luttinger liquid 
regime V -C l/r c . 



/( 3 ) r» £/3y3g(l-a)-2 T 3g(l-a)^ ^ 

At small U the second order contribution always ex- 
ceeds higher-order contributions. Interestingly, however, 
at greater U which are still within the region of validity of 
the perturbation expansion, the third order contribution 
is leading. 

Let us show that the third order contribution domi- 
nates. We will focus on the largest Ui accessible with 
the perturbation expansion: Ui ~ l/r c at g(l — d) > 1 
when Ui is irrelevant; and Ui ~ (VT c ) 1 ~ 9 ( 1 ~ a > /t c at 
g(l — a) < 1 when Ui is relevant. In all cases the elec- 
tron interaction is repulsive, i.e., g < 1. Let us first 
consider the case of g(l — a) > 1. According to Eq. (|23|) . 
a > —1/3. Hence g(l — a) < 4/3. The ratio of the sec- 
ond order contribution to the third order contribution 
scales as (Utc)- 1 (Vt c ) 2 - 9 ^-^ ~ (Vt c ) 2 - 9 ^-^ < 1. 
Thus, the third order contribution dominates indeed. 
The case of g(l — a) < 1 is also easy. In the limit 
of Ui ~ (Vtc) 1 -^ 1 -® /t c we find t c I^ ~ (Vt c ) 2 and 
/( 3 ) V ^> I^ 2 \ Interestingly, the third order current 
may become comparable to the total ac current ~ e 2 V/h 
through the junction at g(l — a) < 1. Note also that 
the exponent 3g(l — a) — 2 in the voltage dependence of 
the dc current (|39|) is negative at g(l — a) < 2/3. This 
corresponds to the dc current increase as the ac voltage 
decreases. 

The above calculation applies for the voltage inter- 
val l/r c > V at g(l - a) > 1 and l/r c > V > 
{UtcY^-^-^/tc at g(l - a) < 1. The left inequal- 
ity is dictated by the applicability of the Luttingcr liquid 
model. The right inequality is determined by the va- 
lidity of the perturbation theory. One cannot calculate 
the rectification current at lower voltages with the above 
perturbative approach. However, it is obvious that 7 = 
at V = 0. From this we conclude that there is a bump 
on the voltage dependence of the rectification current at 
g(l — a) < 2/3: it grows as a function of the voltage at 



IX. RECTIFICATION EFFECT IN THE 
PRESENCE OF THE TIME-REVERSAL 
SYMMETRY 

This case is more complicated than the situation with- 
out the symmetry. The summary of the results is given 
in section IX. C. We will need to take into account many 
different backscattcring operators. Thus, it is important 
to classify them. The classification relies on the matching 
conditions at the junction. 

The matching matrix is given by equation (Ti"2"|) where 
the expressions for a = fi{a),b = /2(a), c = /3(a) in 
the parametrization^ 2 - are written under the equation. 
It is always possible to redefine a in such a way that 
a = fk(a),b = fi(a),c = /„(«), where (k,l,n) is an 
arbitrary transposition of (1, 2, 3). Indeed, the change of 
the variable a — J- 1/a exchanges fi and f2", a — >• — (1 + 
a) exchanges /1 and /3j a — > —a/(l + a) exchanges 
/2 and fz. Any other transposition is a superposition 
of the above three. Since M is an orthogonal matrix, 
at least one of the elements must be negative or zero. 
Otherwise, different rows cannot be orthogonal. Without 
loss of generality we may assume that c < and hence 
a > 0. This can always be achieved by renumbering the 
wires and redefining a. Note that a, 6 > 0. Similarly, 
without loss of generality we may assume that a > b. 
Hence, a < 1. Note that a > 2/3. Thus, the ranges 
of the parameters that we consider are: < a < 1, 
-1/3 < c < 0, < b < 2/3 and 2/3 < a < 1. 

As has already been discussed, we need to deal with 
two types of backscattering operators in the action and 
the current operator: 

(*) Ui - Vexp(in k [<fe(x = 0) - = 0)]); 

(**) Wi - U 2 expiinilfyix = 0) - 4>j(x = 0)] + 
in a [<fi } (x = 0)-<ft(x = 0)]). 
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Note that we do not need to include a contribution of the 
form 7T-2 \<po (x = 0) — 4>~i(x = 0] in the exponent in Wj. 

We will first concentrate on the operators of the 
first type and determine the leading contribution to the 
ratchet current in the absence of the operators of the 
form Wi. Next, we will check what changes after Wi are 
taken into account. We summarize our findings in Sec. 
IX. C. 



A. Contributions from operators Ui 

The scaling dimensions of the operators Ui are simply 

zi = n\g(l - b):z 2 = njg(l - c); z 3 = n%g{\ - a). (40) 

The lowest scaling dimension is z 3 with 713 = 1 . It sets the 
scale of the maximal U in the perturbation theory, U ~ 
(W c ) 1-fl, ( 1_ °) /t c . Strictly speaking we need to know the 
scaling behavior of the operators Wi to determine the 
maximal allowed value of U in the perturbative regime. 
We will see in the next subsection that the operators 
Wi do not change the above expression for the maximal 
U. We thus immediately find the leading second-order 
contribution to the rectification current 

jPO ~ —(UT c ) 2 {VT c ) 29 ^- a) . (41) 

What about the third and higher orders? It turns out 
that the leading contribution is third order and comes 
from the operators U exp(2i[<pQ(x = 0) — ${x = 0)]), 
t/exp(— i[0^,(x = 0) —$(x = 0)]) and Uexp(— i[<f>%{x = 
0) — 4> 3 j(x = 0)]). It scales as 

jO) „ i([/ Tc )3( FTc ) 6 »t 1 -«l- 2 (42) 

To see that (|42|) is the leading contribution, we first 
note that I^ 3 \ Eq. (|4"2"j) . exceeds i' 2 ' at the maximal al- 
lowed U ~ (Vtc) 1 -^ 1 -^ /t c . Indeed, at such U. I (2) t c ~ 
(Vt c ) 2 while I^t c ~ (Vt c ) 1+3 3^] > (Fr c ) 1+3 9 xl /3 > 
(Vt c ) 2 . Next, let us compare the contribution (|42|) with 
other higher order contributions to the current, 

I' ~ Vj[ (t/V^Vf ) ~V]\ ([VT c ] z '-^- a ^ . (43) 
1 1 

If an operator with the scaling dimension zi = ri\g\V — c] 
enters the above expression then (|43|) is smaller than 
(|42j) . Indeed, any contribution to the rectification cur- 
rent that contains such operators satisfies the inequal- 
ity I't c < (VT c ) 1+ ^ a -^ = (y Tc )l+9[l+2a]/[l+"+a 2 ] < 

(V r r c ) 1+39 [ 1_a ] . Thus, only contributions with the op- 
erators Ui ~ U exp(inu[(j)Q(x = 0) — ^(x = 0)]) and 
U[ ~ U exp(in 3 i[(j) 3 :) (x = 0) — <f) 3 (x = 0)]) have to be con- 
sidered. If one of the contributing operators has \nu\ > 1 
or |n 3 ;| > 1 then it is obvious that the contribution 
cannot exceed (|4"2"1) . We thus consider the case with all 



nu,n 3 i = ±1. A nonzero contribution to the current re- 
quires (£j "11 )[#>(* = 0)-(ft(x = 0)] + (£, »3l)[$>(s = 
0) — 0/ ( x = 0)] = 0- F° r a general transcendental a this 
equality is satisfied only if rii/ = 77,3; = 0. Thus, 
we consider contributions in which the operator U\ = 
UiF^F} cxp(iai + i[4>o{x = 0) - <f>}(x = 0)]) enters the 
same number of times as the operator Jj\ and the opera- 
tor U 3 = U 3 F^Ff exp(ia 3 +i[(f> 3 {x = 0)-<ft(x = 0)]) en- 
ters the same number of times as the operator . Hence, 
the phases 011,3 drop out from the final answer. Let us 
now change the voltage sign, perform the transformation 
4>jiq — > —4>i/oi change the order of the Klein factors 
in each term and redefine Fjj —J- FjJ Q , Fjj —> Fj^ . 
Changing the order of the Klein factors introduces com- 
plex conjugate phase factors into Uu and U^.. Hence, we 
can ignore both those phase factors and ctk since nei- 
ther affects the final result. On the other hand, if we 
ignore the phase factors, we discover that the backscat- 
tcring part of the Hamiltonian docs not change under 
our transformation while the current operators change 
their signs. This means, in turn, that the contribution 
to the total current we are calculating changes its sign 
when the bias voltage changes its sign. Hence, it does 
not contribute to the rectification current. Thus, Eq. (|42p 
describes the main contribution to the rectification cur- 
rent. 



B. Contributions from operators Wi 

We now find the leading contribution to the rectifica- 
tion current that contains operators Wj. First of all, let 
us check that renormalized operators Wj(^) ~ Wi(E ~ 
l/r c )(V A r c ) Zi ~ 1 remain small at U < (Vtc) 1 -^ 1 -^ / V c . 
Indeed, Wi(1/t c ) ~ U 2 and z t > 0. Hence, Wi(V)t c < 
{Vtc) 1 -^ 1 -^ < (Vtc) 1 - 2 ^ 3 < {VT c f/ 3 < 1. Thus, all 
renormalized operators remain small and the maximal U 
for which the perturbation theory can be used was found 
correctly in the previous subsection. In what follows we 
focus on the case with U of the order of its maximal al- 
lowed value. The opposite limit U — > is trivial since 
in that limit the second order contribution in U always 
dominates the ratchet current. 

Next, let us check how the operators Wi affect sec- 
ond order contributions to the rectification current. 

(2) 

The corresponding second order contribution ± w ~ 
[Wi{V)fV 2 /rl < (Ut c ) 4 /t c ~ (Vt c ) 4 -^-]/ Tc . Wc 
need to compare it with the current (|42| . ~ 
(W c ) 1+39[1 " ol /Tc- One easily sees that 4 - 4g[l - a] > 
1 + 3g [1 — a] since g[l — a] < 1/3. Thus, the second 
order contribution to the rectification current can be ne- 
glected in comparison with the third order contribution 
even after Wj are taken into account. 

At the same time, in a certain region of parameters 
the dominant higher-order contribution to the rectifi- 
cation current contains an operator Wi. As we will 
see, such contribution contains exactly one operator Wi. 
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We thus begin our analysis by excluding contributions 
which contain three or more Wi's. Indeed, any such 
contribution I 3W < V(Wi/V) 3 ~ (Ut c ) 6 /(V 2 t c ) ~ 
(V^T c ) 4_6£, [ 1 ~ a ]/r c . We need to compare this estimate 
with /( 3 ) - (Vt c ) 1+39 ^-^/t c . One easily sees that 
[4 — 6.g(l — a)] > [1 + 3.g(l — a)] and hence we do not need 
to take into account contributions with three or more 
operators Wj. 

Thus, it remains to consider higher-order contributions 
with one or two operators Wi- Wc first consider contri- 
butions with exactly one operator Wi = W\ exp(ia; + 
J2k s k[4>o( x = 0) — (/>i(% = 0)]), where s 2 can be set 
to zero as discussed above. The contribution also con- 
tains at least two operators Ui. It is sufficient to consider 
the case with all rik = ±1 in the definition of the op- 
erators Ui (*). Indeed, a contribution with an operator 
with I rife I = p > 1 can be increased by substituting it 
with p operators Uk,± = exp(±i[4>^ — <f>o\). Moreover, we 
can assume that for each k only operators Uk,+ or only 
operators Uk,- enter. Let us denote the number of the 
operators of the form exp(±z[</>J — 4>o\) m the expression 
for the contribution under consideration as Ik- Then we 
can estimate the contribution to the rectification current 
as 



j( 3 ') ^ v(U 2 V z ~ 1 )[UV 9{1 ~ a) ~ 1 } h [UV 9[1 ~ b] ~ 1 } h [UV 9{1 ~ c) 

(44) 

where z = g[sza. + si] 2 /(l + a + a 2 ). The same argument 
as in the previous subsection shows that the maximal 
contribution corresponds to I2 = 0. Then for a general a, 
the product of the vertex operators Ui , Wi gives a nonzero 
result after averaging with respect to the quadratic part 
of the action only if l\ = |si| and I3 = | S3 1 . Note that we 
can assume that both si and S3 are nonzero and Si =/= S3. 
Otherwise the operator W; would have the same from as 
one of the operators Ui and the analysis from the previous 
subsection would apply. 

We want to compare the contributions (|4"2"|) and (pFfj) . 

At U ~ (VtcY-^-^/tc we find 



log[/( 3 'V/ (3) ] =7(*i J a3)log(Vr c ); 

7 = 1+ 9 [-5a 2 + | Sl |(l - a 2 ) + ( Sl + s 3 a) 2 }. 
1 + a + a z 



(45) 



If 7 is positive for every choice of s\, S3, si — S3 7^ then 
/(3) is the main contribution to the rectification current. 
If 7 < for a certain choice of Si,S3 then the leading 
contribution comes from J^ 3 \ We thus want to inves- 
tigate at what conditions 7 < and find what choice 
of si and S2 minimizes 7. That choice determines the 
power dependence of the rectification current P3|) on the 
voltage. 

First, let us prove that 7 is minimal, if |si| = 1. Indeed, 
let us compare 7(si = 1,S3 = —1) with 7(^1,^2), where 



\pi I > 1 and pi is arbitrary. One finds 
l{Pi,P2) - 7(1) -1) 



1 + a + a 2 
9 



[(|pi| - 1)(1 - a 2 ) + (p! + ap 3 ) 2 - (1 - a) 2 



~ l + a + a 2 {1 ~ " + a) ~ (1 " a)] > 0l 

(46) 

Thus, we can focus on s± — 1 (the case of si = —1 is 
completely analogous) . 

Next, we prove that a negative 7(1, S3) is minimal at 
S3 = —1. Indeed, a negative 7(1,53) implies that 1 + 
g[l -6a 2 ]/[l + a + a 2 } < 0. Hence, a > (l + v / 41)/10 > 
0.7. Taking into account that a < 1, we find that (si + 
S3<a:) 2 = (1 + s 3 a) 2 is minimal at S3 = — 1. This allows 
one to establish that 7 is minimal at S3 = — 1. 

The remaining task is simple. We just determine at 
what conditions 7(1, —1) is negative, i.e., we need to in- 
vestigate the inequality 



7(1,-1) = 1 + 



g[2 -2a- 5a 2 } 



< 0. 



(47) 



One easily sees that 7 can only be negative, if g > 3/5, 
i.e., for relatively weak repulsive electron interaction. 7 
is negative in the largest interval of a at g —J- 1. In 
-liljhat case, 7 < for a > 3/4. The contributions with 
Wi matter only if 7(1,-1) < 0. They give rise to the 
rectification current of the form (|52|) at negative 7. 

The last question we must address in this subsection 
concerns the role of the contributions to the rectification 
current with two operators of the form (**). An estima- 
tion of such contributions is similar to (l44ll: 



j( 3 ") v(U 2 /V) 2 V z * +Zy [UV 9( - 1 ~ a ^ 1 } h [UV 9 { 1 ~ b ^ 1 } h , 

(48) 

where z XtV = g[s^' v a + s\ ,v ] 2 / (1 + a + a 2 ) are the scaling 
dimensions of the operators W x y ~ U 2 exp(is^' y [4>q{x = 
0) - <t>]{x = 0)] + is x ^[<j>%{x = 0) - <t>j(x = 0)]). 

As above, we first divide I^ 3 ^ by I 1 - 3 * 1 and ask when 
the ratio is greater than 1: 

log[j(3")//(3)] = ^,4,^,4) \og(Vr c ); 



1 + a + a' 



■(-7a 2 + h(l- a 2 ) + (*i + ^«) 2 )- 

(49) 



r—x^y 



If fx < then I^ 3 ' > I^ 3 \ Obviously, a negative 
implies that > 2 — Iga 2 j(\ + a + a 2 ) and hence a > 
(1 + vTL~)/5 > 0.8. Let us now compare fi and 7(1, —1) 
at a > 0.8. One finds 



/i- 7(1,-1) > 1 + 



2 9 



1 + a + a- 



■{a -a 2 -I). (50) 



One easily checks that the above difference is positive at 
a > 0.8. Hence, contributions with one or no operators 
W\ are always more important than contributions with 
two such operators. 
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C. Summary for systems with time-reversal 
symmetry 

We found that the leading contribution to the rectifi- 
cation current depends on the sign of 7(1, —1), Eq. (|47| . 
For strong repulsive interaction (g < 3/5), 7(1,-1) is 
always positive. If 7 > 0, at sufficiently low voltages 
V ~ V* = (t c [/) 1/[1 ~ 9(1 ~ q)1 /t-c the current scales as 



6gc,* 



I r ~ + ° 

At a negative 7 and low voltages V ^ V* 



Ir 



u 4 v- 



^ — 



(51) 



(52) 



The exponents in both voltage dependences are nega- 
tive. This is related to the fact that our calculations are 
only valid in an interval of low voltages, l/r c ^> V > 
V* = (T c U) 1/[1 ~ 9il ~ a)] /t c . Similar to Section VIII the 
I-V curve for the rectification current exhibits a bump 
at V ~ V*. The maximal rectification current at 7 > 



(53) 



At 7 < 



e 2 V 



(54) 



One can easily verify that the maximal current 
()53I54[) exceeds the maximal possible rectification cur- 
rent at the same voltage F in a linear wirei£ ~ 



^max([^T c ] 3 M^T c ]). 



X. CONCLUSIONS 

We have found the rectification current in the absence 
(j39]) and in the presence (|51|52| of the time-reversal sym- 
metry in Y-junctions. In all cases the maximal rectifica- 
tion current is greater than in a linear wire with the same 
interaction strength and bias voltage. In the absence of 
the time-reversal symmetry the rectification current can 
be comparable with the total ac current ~ e 2 V/h for suf- 
ficiently strong interaction strength, i.e., it achieves its 
maximal possible order of magnitude. This reflects the 
fact that both electron interaction and time-reversal sym- 
metry breaking facilitate rectification. For most values of 
parameters the rectification current is a nonmonotonous 
function of the bias voltage. 

Our calculations are valid in the vicinity of various 
fixed points in the low voltage regime. In a general case, 
a junction is controlled by a stable fixed point at low volt- 
ages. For repulsive interaction of spin-polarized particles 



there is only one stable fixed point: three disconnected 
wires^ We found a stronger rectification effect near that 
fixed point than for two weakly connected wires. How- 
ever, the current is low for weakly connected wires. The 
diode effect is much stronger in the vicinity of unstable 
fixed points. Thus, it is important to understand how to 
tune the system close to those fixed points. Some of them 
may be tricky to realize experimentally. Indeed, as dis- 
cussed above, negative elements in the matching matrix 
M imply Andreev reflection and could be obtained in a 
hybrid normal-superconductor structure or in cold atom 



systems 



At the same time, it is straightforward to 



make an "island junction"— , Fig. 1, with positive match- 
ing matrix elements 



M = 




(55) 



For example, one can use three line junctions between 
three quantum Hall systems. The role of impurities is 
played by three constrictions in the junctions. Alterna- 
tively, one can use a single quantum Hall island confined 
between edges AB, CD and EF (Fig. 1). The rectification 
current scales as 



o 12 3 2 



(56) 



at V ~ V* = (Ut c ) < 3 -s)< 1 -9) /t c . Time-reversal sym- 
metry is broken in such setup. Thus, the rectification 
current can be made comparable to the total ac current 
even in the low voltage regime. This is the main result 
of the paper. 

Most of the time we ignored phase factors exp(iai) 
in tunneling operators. As the example of three weakly 
connected wires shows, for special values of the phases 
the rectification effect is suppressed. It has the same 
order of magnitude for other values of the phases. 

We considered the simplest example of a junction: 3 
spin-polarized wires. It would be interesting to general- 
ize our results to a system with spin. This may result in 
a more complicated behavior as the minimal model with 
spin includes 6 channels: two for each wire. Still, we 
expect a similar physics. In particular, a time-reversal- 
invariant system with spin can be obtained from two 
copies of the "island junctions" with opposite spins and 
chiralities. In the absence of the interaction between the 
copies, the problem reduces to a spin-polarized island 
junction. An actual realization, based, e.g., on a topolog- 
ical insulator, must involve interaction between opposite 
spins. We expect that such interaction does not change 
the qualitative picture. 
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